Theory of the isotropic-nematic transition in dispersions of compressible rods 
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We theoretically study the nematic ordering transition of rods that are able to elastically adjust 
their mutually excluded volumes. The model rods, which consist of a hard core surrounded by a 
deformable shell, mimic the structure of polymer-coated, rod-like fd virus particles that have recently 
been the object of experimental study [K. Purdy et al., Phys. Rev. Lett. 94, 057801 (2005)]. We 
find that fluids of such soft rods exhibit an isotropic-nematic phase transition at a density higher 
than that of the corresponding hard-rod system of identical diameter, and that at coexistence the 
order parameter of the nematic phase depends non monotonically on the elastic properties of the 
polymer coating. For binary mixtures of hard and soft rods, the topology of the phase diagram 
turns out to depend sensitively on the elasticity of shell. The lower nematic-nematic critical point, 
discovered in mixtures of bare and polymer-coated fd virus particles, is not reproduced by the theory. 

PACS numbers: 61.30.Cz, 61.30.Vx, 64.70.Md, 61.25.Hq 
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I. INTRODUCTION 

Recently developed methods to control the contour 
length and the effective diameter of elongated virus par- 
ticles (li (21 have opened up the way to systematically 
study the bulk phase behavior of mono- and bidisperse 
rods over a relatively wide range of lengths and widths. 
This is important, because it allows for the experimental 
verification of a vast amount of theoretical work that has 
been done on the isotropic-to- nematic {IN) phase tran- 
sition in dispersions of mutually repellingrods, and in 
particular that on bidisperse mixtures Ufl IE IE IE 111 - 
The agreement of theoretical predictions for binary dis- 
persions with experimental data |E 0, 0, 
has so far not been as impressive as that for monodis- 
perse ones J£j, where the impact of, e.g ., the size [iTj . 
shape 01 , molecular flexibility 0, IIE |^ and Coulomb 
interactions I^E IIE 13 appear to be well understood 



One of the reasons for this state of affairs is probably 
that in binary mixtures there is a coupling between or- 
dering, fractionation, and demixing, leading to a much 
more complex phase behavior 0, IE IIE HE Hll ■ Indeed, 
in bidisperse systems more length scales compete with 
each other, presumably making them more sensitive to 
the effects of flexibility, residual Van der Waals attrac- 
tions, non-additivity or charges [E HE HE 1^ 113 ■ bi- 
nary mixtures are therefore a much more critical gauge 
of the accuracy of theories than monodisperse systems 
are. In some cases, the disagreement between theory and 
experiment is not just quantitative ^E| but even qualita- 
tive. For example, in experiments on aqueous mixtures of 
bare fd virus particles and fd virus particles onto which 
polymeric chains are grafted, Purdy et al. [2| discov- 
ered a nematic-nematic coexistence region that exhibits 



a lower critical point. Extensions of Onsager's classical 
second- virial theory for inflnitely rigid rods incorporating 
a diameter bidispersity do not predict such a lower crit- 
ical point 4, 27, 28, 33J, not even if one allows for non- 
additivity of the interactions between the species [s^ . 
However, these theories do predict a nematic-nematic 
demixing either with an upper critical point or no critical 
endpoint at aU H IH HUHE. 

Theoretically, nematic-nematic demixing transitions 
of bidisperse rods with a lower critical endpoint have 
been predicted, but only for rods with a sufficiently 
large bending fiexibility |^3|; o^' fo'' sufficiently short 
rods for which higher order virials become important 
EE SHE EE Ei El- it appears, however, that some of 
the predictions are quite sensitive to the invoked approx- 
imations, and in our view the issue remains contentious. 
For instance, depending on how precisely higher order 
virials are approximately accounted for, the lower critical 
point appears and disappears in the phase diagram [isj . 
This is reminiscent of the qualitative difference in the pre- 
dicted phase behavior of hard-rod mixtures, depending 
on whether the exact orientational distribution is used or 
a Gaussian approximation to that |EHEI3HE0|- So, 
given the apparently inherent sensitivity to the details 
of the theories, it remains important to explore alterna- 
tive explanations of the observed phase behavior, and in 
particular that of the lower critical end point. 

A factor that has not been considered, and that could 
significantly influence the phase behavior in systems such 
as those studied by Purdy and coworkers, is the flnite 
compressibility of the polymer coating of the rods. The 
effective diameter of these polymer-coated rods is not 
necessarily flxed but could well be a function of the ther- 
modynamic state of the dispersion. So far, the dimen- 
sions of rod-like colloids have been treated as quenched 
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variables, i.e., as invariants of the thermodynamic state 
of the suspension. (An exception to this are micellar 
rods with annealed length distribution ^^.) With the 
polymer-coated rods of Ref. in mind, it seems op- 
portune to explicitly consider the elastic response of this 
coating to the osmotic pressure of the dispersion. Con- 
trary to for instance theoretical work on the crystalliza- 
tion of dendrimers, where entropic interactions of a simi- 
lar nature are modelled by a soft potential , we model 
the effects of compression not at the interaction potential 
level but at the level of the volume exclusion. In effect, 
our theory is that of particles with an annealed diameter, 
and hence an annealed excluded volume. 

Our calculations show that monodisperse fluids of elas- 
tically compressible rods exhibit an isotropic-nematic 
phase transition at a density higher than that of the cor- 
responding incompressible-rod system of equal diameter, 
as expected. How much higher, depends on the stiffness 
of the shell. Weakly impacted upon by the shell stiffness 
is the relative density gap in the coexisting phases. For 
binary mixtures of incompressible and compressible rods, 
the structure of the phase diagram changes dramatically 
with varying elasticity of the coating, showing once more 
the sensitivity of this kind of mixtures to the details of 
the interactions js^l . Still, the topology of the phase di- 
agrams that we calculate does bear some resemblance to 
that of incompressible hard-rod mixtures of unequal di- 
ameter: an upper but no lower nematic-nematic critical 
point is produced by the theory. 

The remainder of this paper is organized as follows. 
In Sec. ^] we first introduce the Onsager-type free energy 
functional, and derive from that the basic Euler-Lagrange 
equations describing the orientational and density distri- 
bution of the elastically compressible rods under condi- 
tions of thermodynamic equilibrium. In Sec. lIIII we study 
the behavior of a monodisperse fluid of such compress- 
ible rods. In order to obtain an analytical solution to the 
model we invoke a Gaussian approximation to the ori- 
entational distribution function. This analytical theory 
we compare with an exact, numerical evaluation and find 
fair agreement. In Sec. lIVI we analyse the bulk phase di- 
agrams of binary mixtures of hard rods and elastically 
compressible ones, but now only numerically keeping in 
mind the sensitivity of the binary phase diagram to ap- 
proximations invoked. A summary and discussion of the 
results are presented in Sec. 



II. DENSITY FUNCTIONAL AND METHOD 

We are concerned with the bulk properties of a fluid 
of cylinders of two different species ct = 1 , 2 of diame- 
ter Dc and equal length L in a macroscopic volume V 
at temperature T and chemical potentials We pre- 
sume the limit D^/L ^ to hold, in which case a sec- 
ond virial theory is believed to be exact The "ef- 
fective" diameter Di of the thin rod is determined by 
the bare hard core of the particle Di = A'^"'''^, whereas 



the diameter of the thick (coated) rods is written as 
1)2(7) = A™''^ + 7AP°' with AP°' twice the thickness of 
the soft polymeric shell. The compression factor 7 G [0, 1] 
parameterizes the deformation of this soft shell due to in- 
teractions with other rods in the system. At infinite di- 
lution we expect zero deformation of the thick rigid rod, 
i.e., 7 = 1, and it is convenient to define the limiting 
diameter ratio d = 1)2(7 = l)/Di = l + /^poi J /s^core _ 

The compression of the grafted polymer layer reduces 
the excluded volume, as we will see explicitly below, at 
the expense of an elastic energy. Introducing the rigidity 
fc, we write the elastic energy of a single rod at a compres- 
sion equal to 7 as A;(7^ -I- 7"^ — 2)/2 in units of thermal 
energy kBT. It is inspired by the theory of ideal (Gaus- 
sian) polymers Isi^, where we presume that the grafted 
polymers are not strongly stretched 0. Hence, we ex- 
pect the numerical value of the rigidity k to be of order 
of the number of chains grafted onto each rod. However, 
for simplicity we presume it to be a free parameter. Note 
that for 7 = 1 the elastic contribution to the free energy 
is zero, as it should. 

The total grand potential P2] of the spatially ho- 
mogeneous suspension is now written as a functional of 
the distribution functions Pcr(u), where u denotes the 
unit vector along the axis of a rod. The distributions 
are normalized such that ^ d\ipa{u) — the num- 
ber density of species a at the imposed chemical poten- 
tial. Within the second virial approximation we write the 
functional as 0,^3 
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^ / d\ip„{u){\n[p„{u)v„] - 1 - I3p.„^ 
+ ^^J dudW E„„, {u; u')p„{u)p„, {W) 

<7<7' 

+ h{^'+^-'-2) J dnp,iu), (1) 

with f3 — (ksT)^^ the inverse temperature, Va- the ther- 
mal volume of the species a, and the excluded volume 
due to hard-core interactions 

E^^, {u, u') ^ L^{D, + D,,)\sm{ai-ccos{u ■ (2) 

where additional 0{LD^) terms are being ignored, in line 
with the needle limit {DajL 0) of interest here. 

The equilibrium conditions on the functional, 
Sn[{p„}]/Spa{u) = and dn[{pa}]/dj = 0, lead to the 
set of nonlinear integral equations. 
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ln[p„{u)i^„]+Y^ J du'E„„,{u,u')p„,{u') 
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fc(7-7 3) J dup2(u) 

2^ J P't(u)p„'(u), 
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to be satisfied by the equilibrium distributions. The 7- 
derivative in the last line involves the 7-dependence of 
D2 through Eq. These equations can be solved, ei- 
ther approximately within a Gaussian approximation or 
numerically on a grid of orientations. Details of the nu- 
merical schemes have been discussed elsewhere j^H • 
Here, in order to find the bulk uniaxially symmetric dis- 
tributions Pa{0i), with 9 = arccos(u-n) the angle between 
the rod unit vector u and the nematic director n, we use 
a nonequidistant 0-grid of Ng = 30 points 9i G [0, 7r/2], 
1 < i < Ng, with 2/3 of them uniformly distributed in 
[0,7r/4]. The resulting distributions can be inserted into 
the functional to obtain the grand potential —pV, with 
p the pressure. Then the complete thermodynamics can 
be inferred, as well as the phase diagram. 



III. MONODISPERSE SOFT RODS 

As demonstrated in Refs. 0)I0|) the IN transition in 
suspensions of rigid rods is a result of a competition be- 
tween the orientational entropy and the packing entropy 
(free volume). In order to estimate the impact of the elas- 
tic term of Eq. (QJ on the IN transition, we first restrict 
our attention in this section to a purely monodisperse 
system of coated rods. Formally, this can be achieved by 
considering the limit /3/ii —00 and pi(u)L^Di ^ in 
the functional and its minimum conditions, and in this 
paragraph we drop the species index "2" for convenience. 

First, instead of numerically solving for the mini- 
mum conditions we adopt a Gaussian Ansatz for the 
one-particle distribution function in the nematic phase, 
with p(u) oc exp(— 00^/2) for 9 < tt/2, and p{u.) oc 
exp(— a(7r — 9)^/2) for 9 > tt/2. Here, a denotes a vari- 
ational parameter that we fix by minimizing the free en- 
ergy. From Eqs. lj2Jl we obtain the following relations 
between the density and the compression in the isotropic 
phase /, where a = 0, and in the nematic phase N, where 
a = 4c^/7r in terms of the dimensionless concentration 
cn defined below: 
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with C = 21n27r"i/2 -|- 3, C/(Ar) = {Tr/4)ni(^N)L^D{l) 
the dimensionless number density of the / (N) phase, 
D(7) = A™'''=-|-7AP°' the effective diameter of the coated 
rods, and £)(!) = A""''" + Ap°K In addition, the dimen- 
sionless pressures p* — {n / A)[3pL'^ D{1) of the isotropic 
and nematic phases can be written as 



p| = c/ 1 -f- C/ 



D{1) 




P*N = 3cAf. 



(5) 



FIG. 1: (a) Equation of state for a fluid of monodisperse 
soft rods {d = i + A""'- / A""''" = 3) in terms of the dimension- 
less bulk pressure p* = {■K/4:)f3pL^ D{1) and the dimensionless 
density c = / A)nL^ D{1) for several values of shell rigidities 
k — 0.1, 1.0, 10. The inset shows the dependence of the com- 
pression factor 7 on density c for the same values of k as in 
the main plot, (b) Bulk phase diagram for the same system 
in c — A: coordinates. The coexistence region (grey area) sepa- 
rates regions of the isotropic (/) and nematic (N) phase, and 
the tie-lines connecting coexisting phases are vertical. The 
dotted lines indicate the corresponding coexistence IN curve, 
calculated with the Gaussian trial function (Eqs. (|KJ). 
The inset shows the dependence of the compression factors 
7/,jv on density c at IN coexistence for the same range of k 
as in the main plot. 



Note that in the isotropic phase (i) the compression 
7/ decreases monotonically with increasing c/, and (ii) 
the elastic contribution to the free energy grows with in- 
creasing density, as one might in fact expect. More inter- 
estingly, in the nematic phase the rigidity and geometric 
parameters of the rods fully determine their compression, 
which is independent of concentration (at the level of the 
Gaussian approximation). Hence, in the nematic phase 
the contribution of the elastic compression can be con- 
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FIG. 2: The same bulk phase diagram as in Fig.^b) in terms 
of the scaled density c-, = {■K/4)nL^ D(y) versus rigidity of 
the polymeric shell k. The right inset shows the dependence 
of the compression factor 7 on the density c~, for the same 
values of k as in the main plot. The left inset shows the ratio 
r = Cj,N/oy,i of the densities of coexisting A*' and / phases 
(solid line) and the nematic order parameter 5* (dashed line) 
as a function of k. The right inset gives the compression 7 of 
the rods in the coexisting phases, as a function of c^. 



sidered as an effect of an uniform bulk field that renor- 
malizes the value of the chemical potential. This is also 
consistent with expression for the pressure of the nematic 
phase which is only indirectly affected by the elasticity 
of the polymer coating of the rods. Overall, the elastic 
term shifts the IN transition to higher densities due to 
reduction of the effective diameter D^j). The results of 
our calculations are shown in Fig. ^ for the system with 
d — 3, and other values of d lead to the similar phase 
diagrams. 

In Fig. nja) we show dependence of the dimensionless 
bulk pressure p* on the dimensionless number density c 
(i.e., the equation of state) for several values of rigidity k 
of the polymeric shell of the rod. The solid lines represent 
direct numerical solutions of Eqs. The results of 

calculations within the Gaussian approximation are quite 
close to those from our numerical calculations, albeit that 
they overestimate the rod densities in the /, N phases at 
coexistence (not shown here for the sake of clarity). 

Imposing conditions of mechanical and chemical equi- 
librium between isotropic and nematic phase, we calcu- 
lated the densities and the compression factors of the 
rods in the / and N phases at coexistence, and present 
these in Fig. ^b) for the realistic range of rigidities of 
k G [0.1,10]. The dotted lines give the results of the 
calculations within the Gaussian approximation, whereas 
the numerical solutions to Eqs. Q are represented by the 
symbols, with the solid lines serving as a guide to the eye. 
Whilst the predictions for the concentration of rods in the 
isotropic phase at coexistence are in almost quantitative 




. pol , . core 
A /A 



FIG. 3: Bulk phase diagrams in terms of the dimension- 
less density ca = (7r/4)nL^ A''"' versus relative thickness of 
the polymeric shell for rods with several rigidities 

k — 0.1, 1.0, 10. The inset shows the dependence of the com- 
pression factors 7 in coexisting phases on ^p°^ / for the 
same values of k as in the main plot. 



agreement, those for the phase gap are not: the Gaussian 
approximation overestimates the phase gap by a factor of 
about two. For rigidities of order fc « 10 the limiting den- 
sities c/ N of the co-existing phases approach the values of 
the corresponding hard rod system. The polymeric shells 
of the rods are then only slightly deformed (7 ~ 1). On 
the other hand, in the limit of small k, the densities at 
coexistence approach those of hard rods with the smaller 
hard-core diameter A'^"'"'^, because in that case 7 — > 0. 
The leveling off occurs only for very small values of k, 
and is not shown in Fig. ^b). The compression factors 
7/_Ar of the rods in the coexisting phases are shown in the 
inset. The values calculated within the Gaussian approx- 
imation are similar to those obtained from our numerical 
analysis. Again, in the isotropic phase the agreement is 
almost quantitative. 

It is clear from Fig. ^ that the dimensionless density 
c at the IN transition increases with a softening of the 
polymer layer, which of course is not all that surprising. 
From a theoretician's point of view, c — (7r/4)nL^Z?(l) 
is indeed the preferable concentration scale, because it 
does not depend explicitly on 7. The question now arises 
whether experimentally the bare diameter D{1) of the 
rods might be determined from the actual concentration 
of particles at which the nematic phase appears |^ i25l | . 
The answer appears to be negative, because the actual 
diameter will not have this value at the point where the 
nematic phase becomes stable. If D{1) were determined 
independently, e.g., from the compressibility at low con- 
centrations, then this would allow one to obtain a value 
for k. 

Theoretically, the structure of Eqs. (jSJ hints at the 
usefulness of the scaled densities Cj — {'K/4:)nL'^D{j), 
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FIG. 4: Bulk phase diagrams of binary mixtures of hard (A^"' = 0) and soft rods d = 3.0 (a), 4.0 (b),4.3 (c), 4.5 (d) in terms of 
the dimensionless pressure p* = (7r/4)/3j3L^A'^°'''^ and the composition x = n2/(ni +n2) for different rigidities: k = 0.1 (sohd), 
k = 1.0 (dashed), k = 5.0 (dotted, only in (c)), = 10 (dot-dashed). In (c,d) A indicate the isotropic I and nematic Ni,N2 
phases of different composition at coexistence, and the A'^A'' remixing point in (c) is marked with *. 



and it is instructive to compare both density represen- 
tations for it helps explain the underlying physics. The 
phase diagram of Fig. predrawn in terms of {cj,k} is 
shown in Fig. |21 It suggests that the density variation 
with k at the IN transition mainly comes from the re- 
duction of the effective diameter D{'-f) the rods. This 
conclusion is supported by the similarity of the values of 
the compression 7/_jv of the polymer shells of the rods in 
the co-existing phases, the dependencies of which on the 
density are presented in the right inset in Fig. |21 

Finally, several other quantities can be useful for qual- 
itative comparisons with experiments. In the left inset in 
Fig. 121 we present the density ratio r = Cj^n/c^j of the 
/ and N phases at coexistence, as well as the nematic 
order parameter S" as a function of the rigidity k. The 
dependence of cn/ci of the unsealed densities on k is 
similar to that of c^^m/c^,i, and is not shown here. The 
small variations of c^^n/c^,i and S with k indicate that 
in experiments it would be quite hard to use these two 



quantities to determine the shell rigidity k. 

Variation of the molecular weight of polymer that 
forms the shells of the rods allows for a modification 
of its thickness Ap°'. In Fig. |3 we show several bulk 
phase diagrams in terms of the dimensionless densities 
CA = (7r/4)ni2A=°'-« and the ratio I^voi j ^core gyg_ 
tems with rigidities k = 0.1, 1, 10. Note that our previous 
definition of the scaled density c — (7r/4)nL^Z3(l) (as in 
Eqs. lOl and Fig.^ would be inconvenient here, as 
it contains explicit dependence on Ap°' . At the smallest 
studied thickness of the polymeric shell /\poi- / /\core — q i 
densities c^jm closely approach correspondent values for 
monodisperse hard rods. As thickness of the shell in- 
creases, the transition densities cajn decrease as one 
would in fact expect. The inset in Fig. O shows the vari- 
ation with the shell thickness of the compression factors 
7/. TV for the same A:'s as in the main plot. It appears that 
the compression of the shells increases for larger values 
of AP°'. 



6 



1 




Y 0.6 



0.4 - 




X 

FIG. 5: Compression factors 7 as a function of the mole frac- 
tion X in the coexisting isotropic (lower curve of each pair) 
and nematic (upper curve of each pair) phases in the bi- 
nary mixture with d = 4.3. Different pairs of curve corre- 
spond from bottom to top to increasing values of the rigidity 
k = 0.1, 1.0,5.0, 10. See also Fig.glc). 



IV. BINARY MIXTURES 

For binary mixtures of soft rods the effective diame- 
ters D„{'y) are different, and one can expects significant 
modifications of the phase behavior in comparison with 
the hard-rod binary fluids. We study the bulk proper- 
ties of mixtures of bare hard and coated soft rods with 
d = 3.0,4.0,4.3,4.5 for various values of k, chosen to 
mimic the experiments of Ref. 2J . Before presenting re- 
sults of our calculations, it is useful to recall the struc- 
tures of the bulk phase diagrams of hard-rod mixtures, 
i.e., of mixtures of incompressible rods with different di- 
ameters. For all values of the parameter d they exhibit 
a low-density I phase, which at some intermediate val- 
ues of densities separates into coexisting I and A'^ phases 
with different composition. In mixtures with d < 4.0 
the high-density region has a single N phase, whereas 
for the systems with d > 4.0 this N phase can demix 
into two different nematics iVi, N2 (depending on the 
composition). The diameter ratio d of the species also 
determines whether this NN phase separation persists 
to (arbitrary) high densities (c? > 4.2), or whether these 
A^i, N2 phases remix back 0,111 El. 

The phase diagrams of the "soft-hard" mixtures were 
determined by solving Eqs. Q under conditions of me- 
chanical and chemical equilibrium. We have verified that 
an artifact resulting from the discretization of the an- 
gular degrees of freedom of the rods, and that produces 
a nematic phase of perfect orientational order, does not 
interfere with our calculations |43 |. 

The results of the calculations are shown in Fig. 01 



where we give the bulk phase diagrams of the mixtures in 
terms of the dimensionless pressure p* = (7r/4)/3pL^A'^°'''^ 
and composition x = n2/{ni + 712) for different rigidities 
k. In this representation the tie lines connecting coexist- 
ing state points are horizontal because they correspond 
to the equal pressure condition. There are several fea- 
tures of the bulk phase diagrams which we would like to 
point out: 

(i) For all diameter ratios d of mixtures of hard and 
very soft rods [k ~ 0.1), the pressure p* at coexistence 
varies almost linearly with composition xm of the ne- 
matic phase. Detailed numerical evaluations show that 
p* depends approximately linearly on the total density of 
the mixture n = ni -f n2, which is reminiscent to the be- 
havior of a monodisperse system. This is to be expected, 
because there is only a weak coupling then to the elastic 
response of the soft component. In our calculations, both 
hard and soft rods have a hard core of equal diameter. 

(ii) For a given value of d, the phase gap at IN coex- 
istence is significantly smaller for systems with smaller 
values of k than that for large k. A large phase gap is 
usually seen as indicative of polydispersity effects. In- 
terestingly, a narrow IN phase gap was observed in the 
experiments with PEG-coated fd- viruses 0, much nar- 
rower than to be expected if both components were in- 
deed incompressible. 

(iii) Although a rigidity of A; = 10 seems high and 
should render the mixture close to that of hard rods of 
unequal diameter, as can inferred from our results on 
monodisperse rods in Sec. lIIII a N1N2 phase separation is 
not observed for d — 4.0. We recall that it does manifest 
itself in rigid hard-rod mixtures of this diameter ratio 

28, 33] . Such a strong sensitivity of the stability of the 
nematic phase on the elasticity of the shell is also seen in 
Fig. ^c), where an increase of the rigidity from /c = 5 to 
k = 10 hardly affects the IN transition curves whereas 
it does prevent the N1N2 remixing at high pressures in 
the fc = 10 case. 

Finally, in Fig. \^ we show the compression factor 
li,Nix) of the polymeric shells of the rods in the coexist- 
ing / and N phases for k — 0.1, 1, 5, 10 for d ~ 4.3, as a 
function of the mole fraction x. Mixtures with other val- 
ues of d and k produce similar curves. Upon an increase 
of the relative concentration x of the soft rods, compres- 
sion of the polymer shells increases monotonically, which 
reflects proportionality of the elastic energy to x. Note 
that 7/ (a;) < 'Jn{x) for all the systems studied. Almost 
linear curves ^i,n{x) for extremely soft rods [k = 0.1) 
again indicate effectively monodisperse behavior. 



V. SUMMARY AND DISCUSSION 

In this paper we have explored the bulk phase dia- 
grams of monodisperse compressible rods and those of 
binary mixtures of incompressible and compressible rods. 
Our motivation for this are the experiments by Purdy 
and coworkers, in which polymer-coated fd virus parti- 
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cles were mixed with bare fd in aqueous suspension 
Not surprisingly, our study shows that a pure system of 
coated rods, modeled here by an elastically responding 
excluded-volume interaction, requires a higher number 
density to obtain a nematic phase than that of incom- 
pressible rods of equal diameter. Nevertheless, the den- 
sity ratio of the rods in the coexisting phases remains 
very close to that of the pure hard-rod system, being 
1.274. The same, in fact, is true for the nematic order 
parameter. This implies that these dimensionless quan- 
tities cannot be used to estimate the elastic modulus k 
of the polymer coating. For known values for the coat 
thickness, A^"', and the bare core diameter, A'^°''^, that 
of the modulus k might be calculated from the observed 
absolute density in either coexisting phase, using, e.g.. 
Fig. n if AP°^ / A'"'"'" = 2. 

The situation is very different when it comes to mix- 
tures of compressible and incompressible rods, which 
we have investigated for relative coat thicknesses in the 
range of values equal to d = 1 + A^"'/^™''" = 3-0, 4.0, 
4.3 and 4.5. We find that the topology of the phase dia- 
gram in the pressure-composition representation changes 
quite dramatically in this small range of d-values, with 
0.1 < fc < 10. Upon an increase of the shell rigidity from 
the lowest to the highest value, we find (i) that the degree 
of fractionation at IN coexistence increases, and (ii) that 
the NN binodal, if present, changes from a dome like to 



chimney like, i.e., the upper critical point moves to infi- 
nite pressures at a critical value of k that depends on d. 
For conditions where there is an upper critical point in 
the phase diagram, we did not find a reentrant NN phase 
separation at higher pressures, i.e., no lower critical NN 
point was found. Therefore, the experimentally observed 
lower critical point cannot be explained by a compress- 
ible polymer coating, at least not within an Onsager type 
of approach. 

As it is now clear that the phase behavior of the rods 
depends sensitively on the presence and properties of a 
polymeric coating, we suggest that a realistic model of 
any experimental system involving coated rod-like parti- 
cles should not only correct for a finite bending rigidity 
and/or finite-length, but also for the elastic softness of 
that coating. This we intend to pursue in the near fu- 
ture. 
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